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A STUDY OF SURFACE PRESSURE FLUCTUATIONS 
IN HYPERSONIC TURBULENT BOUNDARY LAYERS 
By K. R. Raman 

Nielsen Engineering & Research, Inc. 

SUMMARY 

An experimental investigation of the surface pressure fluctuations on 
a flat-plate model at hypersonic speeds of 5.2, 7.4 and 10.4 with an 
attached turbulent boundary layer was carried out at Ames Research Center. 
The statistical data obtained include RMS pressures and nondimensionalized 
energy spectra for all the three Mach numbers. Space-time correlations 
and the convective .velocities were obtained for tests at Mach number of 
7.4 only. High frequency resolution of the pressure field was achieved 
using specially designed small piezoelectric sensors (diameter = 0.5 mm) 
that had a frequency response well above 300 KHz. The results indicate 
the RMS pressures vary from 5 to 25 percent of the mean static pressures. 
The ratio of RMS pressures to dynamic pressures is less than the univer- 
sally accepted subsonic value of 6x10”^. The ratio decreases in value as 
the Mach number or the dynamic pressure is increased. The ratio of RMS 
pressure to wall shear stress for M^ = 7.4 satisfies 1 1 . 

broad-band space-time correlations with different streamwise spatial sepa- 
rations are characterized by a convective velocity U^ nearly equal to 
0.7 of the free-stream velocity. 

INTRODUCTION 

The generation of sound is in general due to vibrations of surfaces 
or fluctuating shear flows in a compressible medium. These vibrations 
produce a fluctuating force in the medium and cause density and pressure 
fluctuations that radiate energy outward. If the intensity and frequency 
of these pressure fluctuations fall in the audible region we perceive them 
as sound. Knowledge of the pressure fluctuations on the surface within a 
turbulent boundary layer at various Mach numbers is essential in many 
problems. Although considerable work has been performed to investigate 
such fluctuations at subsonic and supersonic speeds, little information is 



availalbe at hypersonic flows. At subsonic and supersonic speeds and at 
large Reynolds nxoinbers, a major portion of the exterior surface of flight 
vehicles is subjected to pressure fluctuations due to attached and/or 
separated turbulent flow. The intensity of these fluctuating pressures and 
their frequency components 300 KHz expected at hypersonic free-stream 
velocities) could significantly influence the dynamic load analysis of the 
structure and lead to fatigue of the vehicle. 

The fatigue loads caused by the fluctuating acoustic pressures may 
determine design requirements for the skin of hypersonic re-entry vehicles, 
cruise vehicles, or logistic resupply space vehicles. The vibrations and 
noise transmitted through the structure can affect adversely the components 
in guidance instruments, decrease flight crew efficiency and affect the 
protective thermal shield of space shuttle vehicles. For structural 
dynamics and analysis-, the designer is interested in being able to predict 
pressure data related to the aerodynamic parameters of the turbulent exter- 
nal flow over the surface of the vehicle at hypersonic speeds. These 
fluctuating pressure data enable one to assess the intensity of a turbulence- 
generated sound field that radiates its energy from sources in the boundary- 
layer region. For a successful suppression of the unwanted noise thus 
generated, it would be necessary to understand the source characteristics. 

The noise produced by the fluctuating pressure field in a turbulent 
boundary layer can be classified into "far" and "near" fields. The turbu- 
lent fluctuations will (i) radiate a sound field into the external flow, the 
"far" field and (ii) create pressure and shear fluctuations adjacent to 
the surface, the "near" field. If the bounding surface is flexible, 
coupling exists between the wall motion and the external flow field. 

A literature survey on the fluctuating pressure field in turbulent 
boundary layers reveals theoretical treatments of the problem by Philips 
(refs. 1 and 2), Curie (ref. 3), Kraichnan (refs, 4-6), Lilley and Hodgson 
(ref. 7), and others. In addition, several experimental investigations 
concerning incompressible boundary layers in wind tunnels have been carried 
out by Harrison (ref. 8), Willmarth (refs. 9-14), Bull (refs. 15-17), and 
Serafini (ref. 18). Laufer (ref. 19 and 20) considered the supersonic far 
field problem. Chen (ref. 21), Kistler and Chen (ref. 22), Lilley (ref. 23), 
Speaker and Ailman (ref. 24), Murphy et al. (ref. 25) Maestrello (ref. 26), 
and Coe and Chyu- (ref. 27) carried out experiments in the turbulent boundary 
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layer at supersonic speeds up to Mach number 3.5. The pressure field 
investigation of separated flow in the boundary layer due to a step or 
interacting oscillating incident shock wave has been performed by Speaker 
and Ailman (ref. 24) and Coe (refs. 27 and 28). In addition to all the 
above literature, there are some free flight data obtained from a glider 
flight, and from powered airplanes (refs. 29-35) . 

All the above literature indicates that the statistical pressure 
information for the incompressible, attached turbulent boundary layer is 
fairly complete. At transonic and supersonic speeds, measurements of the 
near field pressure intensity have been carried out; in addition there 
exists some spectral and correlation information for the attached super- 
sonic flows. From these limited experimental results at supersonic speeds, 
it can be concluded that the data are consistent with the findings at sub- 
sonic range and that the information on the pressure field is well in hand. 
Beyond the supersonic range there is a scarcity of data and no direct 
pressure fluctuation measurements have so far been reported in the liter- 
ature. Wallace (ref. 36) and Harvey et al. (ref. 37) have given some 
pressure intensity values at hypersonic velocities. These were obtained 
indirectly by calculations from temperature and density fluctuation data. 
These measurements were made with the use of an electron beam technique at 
hypersonic speed. The present investigation of near field pressure fluc- 
tuations at hypersonic speeds is one of the first studies that provide 
direct measurements. 

The experimental research reported here is concerned with the problem 
of a fluctuating pressure field within an attached turbulent boundary layer 
over a flat-plate model when the free-stream velocities are in the hyper- 
sonic range (5 £ M^ ^ 10.4) . The tests were carried out in the Ames 
Research Center 3 . 5 Foot Hypersonic Wind Tunnel, and the range of Reynolds 
numbers in the experiments was from 0.9 to 15 million. 

A major part of the work reported here was carried out at Ames 
Research Center, Moffett Field, CA by the author as an employee of NASA 
prior to October 1971. The data analysis and preparation of the report 
were carried out at Nielsen Engineering & Research, Inc., Mountain View, CA 
on Contract No. NAS2-6838 from Ames Research Center. 
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LIST OF SYMBOLS 

constant related to ratio of RMS pressure to dynamic pressure 
see equation (13) 

local skin-friction coefficient 

diameter of sensor, meters 

frequency, Hz 

low frequency roll-off setting, Hz 
high frequency roll-off setting, Hz 

gravitational acceleration, 980 cm/sec^ 

nondimensional power spectral density function of 
p(t) = p^ (f ) U^/q^6* , see equation (10) 

transfer function of sensor in frequency domain 

free-stream Mach number 

shock wave speed 

static pressure, N/m^ 

energy or power spectral density function of p (t) , (N/m^ ) ^/HZ 

00 

mean square value of pressure fluctuations, = jp^(x,f)df, 
(N/m^) o 

roo t-mean- square (RMS) value of pressure fluctuations, 

~\J p^ (x) ,N/m^ 

reservoir pressure, N/m^ 

pressure step obtained in the shock tube calibration apparatus, 
N/m^ 

free-stream dynamic pressure, N/m^ 

percent variation in Reynolds number, see equation (6) 

Reynolds nxmiber 

normalized auto-correclation function, R..(x,x) = 
11 

p^(x,t)p^(x,t + t) p| 
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normalized cross-correlation function, R. .(x,t) = 

ij 


(x , t) Pj (x, t + T ) /piPj 


normalized space-time correlation function R^^ (x , 5 , x ) = 
p^(x,t)pj(x + |,t + t )/^ p^(x)pj(x + t) 


reservoir . temperature (°K) 


free-stream velocity, meters/sec 

broad band convective velocity, m/sec, given by equation (20) 

broad band convective velocity, m/sec, given by equation (22) 

vector representation of coordinate system, (x^,X 2 ,Xj) 
distance 'from leading edge to instrument plug, meters 
ratio of specific heats 

standard boundary-layer parameters, meters 

vector representation of separation distance between sensors 

fixed separation distance between sensors ( =9.5 mm) 

time delay, sec 

local wall shear stress, N/m^ 

circular frequency, ( = 2iTf ) , radians/sec 


A,B. . .H 
on p 

X , 6 , 6 * and 
6 on Re 

j 


Subscripts and Superscripts 

represents measured pressure from sensor A,B or . . .H, 
respectively 

Reynolds number based on x, 6, 6* or 6, respectively 

assumes any of the sensors; namely, i = A...H and j = A...H 
bar represents time average 
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WIND-TUNNEL FACILITY 


The tests for this investigation were carried out in the Ames 3.5-Foot 
Hypersonic Wind Tunnel. Schematic drawings of the facility are given in 
figure 1(a) and (b) . Air, free of dust and water vapor, from high pressure 
tanks passes through a heated alumina pebble bed chamber, an interchangeable 
nozzle section, a test section and a diffuser portion of the tunnel prior 
to reaching large evacuated spheres. From these' spheres the air is returned 
to the high pressure tanks. The air is preheated by the pebble heat ex- 
changer in order to prevent liquefaction in the test section. 

Because of the large frictional losses due to large velocity gradients 
at the interior wall surfaces, small amounts of helium are injected tangen- 
tially to the interior surface at the throat of the nozzle section. This 
reduces the power required for the wind-tunnel operation and increases the 
useful core cross section of uniform test stream and substantially maintains 
the tunnel walls at a low temperature. After each test run the air is 
cleaned to remove dust and moisture, and the helium is separated by a 
cryogenic separation process. The air is then returned to the.ihigh pressure 
storage tanks. 

The wind tunnel is capable of maintaining a continuous steady flow in 
the test section for any of the preselected nozzle configurations at hyper- 
sonic speeds for a maximum duration of 5 minutes. For each preselected 
test stream Mach number of 5.2, 7.4 or 10.4, the appropriate nozzle confi- 
guration must be installed between the pebble bed heater and the test 
section. The wind-tunnel temperature and pressures were monitored both in 
the reservoir (in the pebble bed itself) and at several locations in the 
test section. 

The tunnel conditions for each of the tests are given in Table I. In 
general the tests spanned the following conditions : 

5.2 < Mach number, M <10.4 

C» 

1/2 _< duration of test run, minutes ^ 1 . 

6.4x10^ £ reservoir pressure, N/m^ ^ 1. 10x10^ 

700 < reservoir temperature, °K < 1200 
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6.8x10^ < dynamic pressure, N/m^ _< e-Gxlo"* 

9x10® < Reynolds number, Re^ 1.45x10^ 

The model was injected into the test stream once the preselected 
conditions were reached in the wind tunnel. The model injector mechanism, 
shown in figure 2, is operated hydraulically. The injector mechanism 
permitted retention of the test model in the stream for a minute or less, 
the time required for data acquisition on a high speed tape recorder. When 
the quick insert mechanism is activated, the side wall of the tunnel slides 
out. The model mounted on the quick insert arm is inserted into the middle 
of the test section at the predetermined angle of attack and the side wall 
closes to prevent the hot test stream from flowing into the quick insert 
chamber. After the data are taken, the model is retracted and the exact 
reverse of the above described process takes place. 

In cases of emergency, the model can be retracted quickly by a built-in 
alternative system combining a spring and explosive charge cartridge. At a 
Mach number of 5.2 and high dynamic pressures, this alternative system was 
used. 

A 

V. V 

TEST MODEL 

The model used in these tests, figure 3, was a 1.9-cm thick aluminum 
plate, 120 cm long and 45 cm wide. The fine dust from the disintegrating 
alumina pebbles in the heater impinged on the leading edge . The wind- 
tunnel dust particles were not a major problem within the boundary layer; 

A hard Inconel leading edge with a nose radius of 0.05 mm^. and internal 
water cooling to withstand the high stagnation temperature was used in the 
tests. The entire surface was honed and lapped to provide a highly polished 
smooth surface. The surface was anodized in order to prevent rough irre- 
gular surface oxidization due to moisture from ambient air. 

The plate was shock mounted on the injector mechanism to prevent any 
gross extraneous wind-tunnel shell, or injector, unit vibrations from 
affecting the measurements. The photographs of the test model in figures 3 
and 4 show front, side and back views of the model with all the items 
described above identified. 
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The plate was provided with two holes of 5.7 cm diameter along the 
centerline, about 73 cm and 84 cm downstream of the leading edge to accept 
the instrument plugs. Underneath the smooth test surface was a completely 
enclosed chamber for the instruments to protect them from the heated air 
and high stream velocities (fig. 4) . The electrical wires from all the 
instruments passed through a double-walled, water-cooled conduit from 
the instrument box and emerged from the injector unit described in the 
wind-tunnel section. 

In order to protect the flow on the smooth side of the plate from any 
spill-over from the underside flow conditions, rigid skirts were provided 
along the entire length of the plate on each side, figure 3. This also 
insured a two-dimensional mean flow condition at the instrument plug. 

INSTRUMENT PLUG 

The instrumentation was concentrated in a plug mounted in the flat 
plate. The ■ instrument plug is a circular 5.7 cm diameter and 1.9 cm thick 
disk which could be installed flush with the plate into either hole along 
the longitudinal axis of the test model. For the present tests the instru- 
ment was installed 73 cm from the leading edge. A noninstrument insert was 
installed at the second hole, 83 cm from the leading edge. The disk with 
various instruments was shock mounted using o' rings as shown in figure 5 
in order to eliminate mechanical vibrations and acceleration effects. In 
the photograph of the fully instrumented plug the sensors are labeled from 
A to H and the two steady-state static pressure sensors are labeled I and 
II. The sensors are^ located 9.5 mm apart center to center in both longi- 
tudinal and lateral directions. 

Each of these sensors was carefully mounted under a bino.cular micro- 
scope, and the annular gap (less than 0.02 mm) between the sensor and the 
instrument plug was filled with silicone rubber to avoid any presence of 
Helmholtz resonance and the accompanying discrete frequency sound. The 
rubber could withstand high' temperatures without appreciable swelling or 
contraction. This precaution was essential in several respects: it 

avoided any roughness causing local stream perturbations and it eliminated 
local effects of Helmholtz resonance and acted as a shock mounting for the 
sensor . 
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The back side of the instrument plug is shown in figure 6. Various 
sensors used for pressure and temperature monitoring and the source- 
follower amplifiers are pointed out. The steady-static pressure reference 
was obtained from the two static pressure orifices by connecting each 
orifice to a coil of 3.5 meters long hypodermic tubing with a 0.5 mm inner 
diameter. The attenuation of the fluctuating pressure components in long, 
small diameter tubes is well known (refs. 38 and 39) and the pressures so 
obtained were fed to a common manifold located on the back side of the 
instrument plug. The manifold supplied the required reference pressure 
for all the diaphragm sensors. The steady-static pressures (prior to 
common manifold) were recorded on tape during the tests. The reason for 
two separate steady-static pressures venting to the manifold was to avoid 
loss of data or sensor due to unexpected clogging of the orifice by dust 
particles. 

■ The thermocouples used were made by thin-film technique and the thermo- 
couple elements were chromel-constantum. Because there was very little 
mass, the thermocouples were good for high frequency measurements. The 
thermocouples T^ to T^ shown in figure 5 were carefully mounted on a Kel-F 
heat insulator to avoid the instrument plug acting as a heat sink during 
the data acquisition period. One temperature sensor (T^ in figure 5) was 
mounted directly on the plate to monitor the temperature rise during the 
tests. 

The sensors mounted on the instrvunent plug are listed below: 

(1) 4 crystal sensors (1 KHz £ freq ^ 500 KHz) 

(2) 4 diaphragm sensors (0.3 KHz ^ freq < 80 KHz) 

(3) 4 thin-film thermocouples mounted on Kel-F heat insulators 

(4) 1 thin-film thermocouple to sense the plate temperature during 

tests 

(5) 2 static pressure orifices (steady-state pressures) 

(6) 1 accelerometer 

(7) 1 reference manifold 
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INSTRUMENTATION 


In these tests several sensors were used with their associated elec- 
tronic amplifiers to obtain signal levels sufficient for data acquisition 
on an analog tape recorder system. The parameters measured were the surface 
static pressure time history at several locations, the acceleration of the 
test model, the surface temperature and the static and total pressure and 
total temperature in the test section. The data acquisition/retrieval 
system, the amplifiers, and various pressure, temperature and accelerometer 
sensors used will be discussed next. 

Data Acquisition/Retrieval System 

A specially designed 20-channel analog magnetic tape record/reproduce 
system made by Ampex (Model 1800H) was used with a recording tape speed of 
3.048 meters/second. The recorder has nine channels of wide band direct 
(400 Hz - 1.5 MHz) record-reproduce and nine channels of FM (DC - 40 KHz) 
record/reproduce modes of operational capability. In addition two channels 
were used for time-code and tape control. An edge track voice channel was 
also available. All necessary additional electronics for reproduce phase 
at two additional speeds of 0.381 and 0.048 meters/second were provided 
also . 

A one-megahertz tape speed reference oscillator signal was recorded 
on a direct data track and was useful in assessing any tape stretch problem 
and in determining dynamic skewness between two adjacent channels. This 
also provided a means of obtaining information on the cross-talk between 
channels. The record and reproduce operations are bi-directional at all 
the selected speeds. The tape recorder can shuttle between any two pre- 
selected start and stop time code information off the recorded data, which 
is useful in the data analysis phase. 

The Time Code Generator (Model 6190) , Time Code Translator (Model 5220) 
and Tape Search and Control (Model 5224) were all made by Astrodata. These 
units were integrated into the tape recorder system as a single unit. 

Shock Tube Calibration Apparatus 

A shook tube was built for static and dynamic calibration of pressure 
sensors and is shown in figure 7. When the aluminum foil diaphragm is 
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ruptured, the atmospheric pressure generates a one-dimensional shock wave 

that travels down the 180 cm long, 5.7 cm diameter low pressure chamber. 

The shock speed is obtained by determining the lapse time between two 

sensors located 5.1 cm apart on the side wall of the low pressure chamber. 

The instrument plug is mounted at the end of the tube. Once the shock 

speed, and the ambient initial test conditions are known, the input 

pressure pulse experienced by the instruments can be obtained from normal 

shock relations (ref. 40) . The shock tube, being portable, was easy to 

use at the wind-tunnel location for periodic calibration checks on the 

sensors. Further, the low pressure chamber could be brought against the 

test model and the suction in the chamber made a perfect seal between the 

shock tube ■ and the model prior to diaphragm rupture. Thus, the entire 

bank of sensors could be checked without seriously disrupting the tests. 

The low pressure chamber was evacuated to 280 mm Hg before the 2. 54x10“ ^cm 

thick aluminum foil failed. In all the calibration tests, the shock speed 

did not exceed M = 3.0. 

s 

The static calibration of the diaphragm sensors was also accomplished 
with the shock tube apparatus. The low pressure chamber of the shock tube 
was evacuated in several steps and held at each pressure for the static 
sensor response to be recorded. For this calibration, the backside of the 
diaphragm sensors were vented to ambient atmospheric pressure while the 
manifold, venting to various sensors, was closed. The output of each 
diaphragm sensor was directly read (without the amplifiers in the circuit) 
through the use of a digital millivolt meter at each low pressure setting. 
This procedure was necessary since the Princeton amplifiers have a high 
pass filter setting of 10 Hz (the lowest operational frequency limit) and 
cannot register any DC voltages. 

Pressure Sensors 

A study of the statistical pressure characteristics of boundary- layer 
turbulence through spectral and correlation (space-time) considerations 
demands a very small sensor with a flat frequency response from low audio 
frequency to ultrasonic frequency. The subsonic and supersonic turbulent 
boundary-layer literature indicates that at hypersonic velocities, pressure 
fluctuations with acoustic energy of interest could be expected up to 
300 KHz (refs. 21 and 22). Since no commercially available transducer with 
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the required frequency response characteristics could be found, a piezo- 
electric sensor was designed and developed for this program (ref. 41) . A 
schematic drawing of the sensor is shown in figure 8. 

The basic element of the sensor is a Lead-Zirconate-Titanate crystal. 
The unique design feature of the crystal sensor is the shock mounting of 
the sensing element. This is essential to avoid a signal being generated 
by mechanical vibration of the plate when the plate is inserted into the 
tunnel. Subsequent tests of the sensor mounted in the plate indicated 
negligible response even when hitting the plate with a mallet. 

The piezoelectric material is ideally suited as an active transducer 
element. The desired qualities are (i) its high sensitivity and linearity 
to applied stresses, (ii) its very high natural frequency characteristics, 
(iii) its wide dynamic range of operation, (iv) its adaptability to sub- 
miniaturization, and (v) its ready availability in disk or cylinder forms. 

In addition, present-day electronics make it possible to use the metal 
oxide semiconductor (MOS) and/or field effect transistor (FET) circuitry 
for the matching preamplifier (Emitter-follower or Cathode-follower circuit) 
with high input impedance and low output impedance. The use of these FET 
amplifiers directly connected to the sensor eliminated the usual parallel 
cable capacitance problem (due to long leads) that degrades the sensor's 
output. This approach also avoided the extraneous noise source generated 
by the motion of the connecting cable. The approach of directly connecting 
a small impedance matching electronics circuit eliminated the need for a 
special low noise connecting cable between the sensor and the FET amplifier. 
A commercially available integrated circuit source-follower amplifier (PCS 
Piezotronics Model 402A11) was connected directly to the sensor in the 
tests (see f ig . 6 ) . 

The choice of Titanium backing material for the piezoelectric sensing 
element provided almost perfect acoustic impedance matching. This is 
essential to minimize the interface interference effects. Further, the 
Titanium rod used in the sensor construction was brought to a point. This 
sharp point provided proper focusing of the incoming acoustic wave to 
progress along the backing material and dissipate its energy at the tip 
(ref. 42). Thus, the usual ringing phenomenon observed in sensors of this 
type was completely eliminated from the sensor used in these tests. 


12 



In addition to the crystal pressure sensors, commercially available 
diaphragm-type Kulite sensors (CPL-070-4) with an overall diameter of 1.78 
mm and a useful frequency range from zero to about 80 KHz were used in 
these tests. These sensors in combination with the crystal sensors made 
it possible to obtain data from very low to very high frequencies. The 
Kulite sensor consists of a stretched thin diaphragm on which four piezo- 
resistive strain gages are vacuum deposited to form a Wheatstone bridge. 

The presence of all the elements on the diaphragm eliminates temperature 
compensation problems. The upper frequency range of the sensor is directly 
dependent on the resonant frequency of the stretched diaphragm. 

The calibration curves for both Kulite and crystal sensors are given 
in figures 9 (a) , (b) and (c) . The small shock tube (discussed previously) 
was used for obtaining both static and dynamic calibration data on these 
sensors-. The solid symbols in these figures were obtained from dynamic 
tests and were usually lower than the static data by about 5 percent.* 
Similar results have been observed by Kendall (ref. 43). The crystal sensor 
calibration data were obtained by dynamic tests only. The calibration data 
of all the sensors used in the tests are given in figure 10, in terms of 
accepted decibel terminology with reference sensitivity given as IV/ybar. 

In addition, the time-rise response of less than one microsecond of the 
crystal sensor was determined by its response to a sudden pressure step 
generated by the nose shock off a bullet (Super 225 Winchester) in flight 
at Mach 3 that swept across the face of the transducer (fig. 11) . From 
the time history of the input and the measured output of the sensor response 
one can determine through Fourier transform methods the transfer function 
for the sensor. The complex transfer function is given by 

H(u) = Fourier transform of sensor response/Fourier 
transform of input pressure step 

The frequency response of both amplitude and phase thus obtained from 
H(w) is given in figure 12 (ref. 44). Shadowgraph pictures of the flight 
of the bullet as the nose shock wave swept the sensor surface were obtained 


*The difference between static and dynamic calibration data depends on 
various parameters, none of which have been systematically established 
to date. It is fortunate that in the present case the dynamic data fell 
systematically lower than the static data; and there is only a small 
percent error incurred in using static calibration data. 
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in order to determine the pressure step, Ap imposed on the sensor. From 
these pictures the incident shock angle and the angle of reflection over 
the sensor's flat surface were determined. Using the known forward velo- 
city of the bullet (equal to shock speed) , the ambient atmospheric condi- 
tion ahead of shock wave, the measured incident and reflected angles, the 
step pressure, Ap imposed on the sensor, was calculated. In the above a 
two-dimensional flow condition was assumed (curvature of the spreading 
shock cone was ignored) and tables from reference 40 were used. The sensi- 
tivity of the sensor obtained in this fashion was identical with the ones 
obtained from calibration with the shock tube. 

To avoid prolonged exposure of the sensitive instruments to the heated 
air stream in the test section, the model was injected into the stream for 
a short time and then retracted (fig. 13) . During the model's stay in the 
tunnel, data were recorded on magnetic tape at 3.048 meters/second tape 
speed. The surface wall temperatures of the model were measured at the 
instrument plug during the model's stay in the tunnel; and data for 30 
seconds prior to model retraction from the stream are given in figure 14. 
The temperature rose from the ambient conditions by less than 10 percent 
of the breakdown Curie temperature of 360°C for PZT-5 material. The 
reservoir temperatures for these runs at the two Mach numbers given in the 
plots were around 1100°C. The temperature drop immediately after the 
model was retracted into the quick insert chamber is noticeable in the 
figure . 


Amplifiers 

Three different types of amplifiers were used in these tests. These 
amplifiers were selected for their compatibility with the sensor and the 
frequency range of interest of the parameter being measured. 

The piezoelectric transducers, being self -generating charge devices, 
exhibit high source impedance to the signal conditioning electronics and 
are essentially all capacitance. Thus, the amplitude and phase characte- 
ristics will be direct functions of the load on the system. In order to 
utilize the excellent undamped high frequency response of the piezoelectric 
sensors and avoid adverse effects on their phase characteristics, the 
associated broad-band amplifier should possess matching high input impe- 
dance and low electrical (or Johnson) noise. The interconnecting cable 
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lengths between the sensor and the matching amplifier should also be kept 
to a minimum in order to avoid degradation of the low frequency characte- 
ristics of the sensor. 

The choice for the tests was a small integrated circuit source-: 
follower amplifier (accomplishing the same functions as a cathode-follower 
or an emitter-follower) -made by PCB Piezotronics , Inc. (Model 402A11) . 

The input resistance and capacitance were 10^ ohms and 10 pf respectively. 
The impedance matching was accomplished without any polarity change, and a 
high-meg resistor drained off any DC components such as those generated 
by temperature. This noninverting amplifier had a voltage gain of unity. 
Thus we needed additional signal-conditioning low noise amplifiers. These 
FET-MOS amplifiers were 6 mm in diameter _ and 2.5 cm long. The source- 
follower amplifier installed in the instrument plug is shown in figure 6. 

The Princeton Applied Research low noise, lock-in amplifiers (PAR 
Model CR-4) powered by rechargeable nickel-cadmium battery packs* were 
selected. These amplifiers were provided with selectable panel switches 
for high and low frequency roll-off settings. In these tests the low 
frequency roll-off was set at 300 Hz and high frequency roll-off at 300 KHz. 
The schematic circuitry for the pressure sensors (Kulite and crystal) is 
given in figures 15(a) and (b) . In order to obtain a low signal-to-noise 
ratio, the source impedance was adjusted by the addition of 2500 ohms in 
series with sensor impedance. Typical amplifier noise figure contours for 
the PAR CR-4 amplifiers are given in figure 16, and the selected operational 
range for the tests is shown in the hatched band in the figure. 

In addition to these amplifiers several wide band differential DC 
amplifiers (Astrodata Model 886) were used for signal conditioning of 
sensors monitoring steady-state quantities like mean static pressure, 
surface temperatures, dynamic pressure, plate location in the test section, 
etc. The measured frequency response of the amplifiers (both PAR CR-4 and' 
Model 886) is given in figure 17. The observed spread between amplifiers 
used in these tests falls within the band shown in figure 17. A mean curve 
is used for the data analysis. The PAR CR-4 amplifier data satisfy 


* 

Battery operation is preferable to rectified 110 VAC power supply. There 
is always some residual ripple left over which cannot be eliminated and 
shows up as noise. 
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Net Relative Gain* = 

[IH- + (f/£„)»] 

where f,f^„ and f,,„ respectively denote the frequencies at which cali- 
J-jr Hr 

bration data were taken, selected low frequency roll-off and the high 
frequency roll-off setting. 


Accelerometer 

An Endevco Accelerometer Model 2225 was mounted on the test model to 
monitor the mechanical vibrations during the pressure fluctuation measure- 
ments. A mounting stud (Endevco #29863) was used to avoid electrical 
ground loop problems. The sensitivity of the accelerometer and amplifier 
combination was adjusted to give 0.707 vrms/1 g. The accelerometer 
traces for a few test runs are given in figure 18. They indicate that 
during initial insertion of the model into the test stream the acceleration 
is less than 1.5 g's and after that time the shock mounting successfully 
isolates the model from the quick insert mechanism. A similar phenomenon 
occurs during the retraction phase. 

Temperature Sensor 

The thermocouple sensors used were manufactured by Heat Technology 
Laboratory (Model TCS-102, with chromel-constantan elements). The sensor 
is 0.8 mm in diameter and about 4 cm long and can be used in an environment 
where the steady-state operating temperature is 1250°K. The time rise 
given by the manufacturer is around 1 psec . 

The TCS Series thermocouple probes consists of a coaxial center lead 
wire (first thermocouple element) housed in a thick wall tube (second 
thermocouple element) . The two metals are insulated from each other by 
a special aluminum oxide-type insulation which is effective to about 


*Net Relative Gain = (amplifier output/amplifier input) x (1/gain setting) 
at the selected signal generator frequency. 

**Source-follower amplifiers similar to the ones used for the crystal sensor 
were used with the accelerometer. In this case an Astrodata Amplifier 
886 was. used since it has an operational range from DC to 20 KHz. 
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1400 °K for steady-state temperatures (and substantially higher for tran- 
sient temperatures) . The thermocouple junction is formed at the polished 
end of the coaxial assembly by a vacuum deposited metallic plating (see 
fig. 19) . Lead wires were welded to the coaxial probe and embedded in a 
refractory insulating cement^. 

For measuring the local stream temperatures the probe was provided 
with a Kel-F (low thermal expansion and low coefficient of conductivity 
material) shielding. Without this shielding the thermocouple sensor would 
be in direct contact with the model and the effects of the model acting 
as a heat sink would be felt in the measurements. 

Correlator 

A Honeywell Model 9410-002 Time Delay Correlator was used in 
obtaining auto- and cross-correlation functions from the recorded random 
pressure time histories of sensors A to H. The Correlator was tested by 
using a sinusoidal signal up to 500 KHz.* 

The Correlator accepts three input time histories simultaneously from 
the tape recorder, say 1, 2, 3, and provides the auto-correlations Rjj(t), 

cross-correlations Rj 2 (t), R 2 3 (t), R 3 j(t). If 
the sensors are separated by fixed distances, the three space-time corre- 
lations are obtained directly when the output of the correlator is recorded, 
namely, R^^(x^,x), R^^(x^,x) and R^j(x^,x). 

The output of the correlator can be recorded with an XY plotter or 
as an oscilloscope trace captured on a polaroid. picture. Both the above 
methods were used for recording the analyzed data. 

Spectrum Analyzer 

For analyzing signals with random contents an estimate of the long 
term power spectrum (or power spectral density characteristics) is usually 
desired. For spectral analysis of the pressure data a Real Time Analyzer 
Model SD301B (Spectral Dynamics Corp.) was used. With its compatible 
Ensemble Averager Model SD302 an estimate of the average spectrum characte- 
ristics of a random signal can be performed up to a frequency of 20 KHz. 

*Even though the manufacturers claim frequency response from 50 Hz to 
250 KHz, the calibration tests carried out prior to our data analysis 
indicated this unit to be good far beyond 250 KHz. 
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The confidence level of the input signal is improved directly in proportion 
to the number of ensemble averages that can be accomplished. With an 
SD302 unit the averaging is accomplished in the frequency domain by 
examining an ensemble of spectral functions, and then obtaining the 
average sum over a predetermined length of time. 

The analysis frequency range of the combination of the SD301B and 
SD302 units can be extended up to 1 megahertz by adding a spectrum trans- 
lator, SD303, and supplying an associated reference signal from a high 
quality Synthesizer (Hewlett Packard 3320A-002) . 

The power spectral plots were obtained with linear horizontal fre- 
quency scale and log vertical spectral density amplitude scale. The output 
is recorded on a Polaroid picture from the oscilloscope display or used in 
an XY plotter . 


Root-Mean-Square Meter 

The root-mean- square values of the recorded pressure fluctuations 
were obtained for each of the eight pressure sensors with a Ballentine 
True RMS meter (Model 320A) which is good for frequencies from 5 Hz to 
4 MHz, and has an accuracy of + 1 percent from the true RMS value. 

TEST PROCEDURE 


While the wind tunnel was being prepared for the test run with the 
preselected reservoir pressure and temperature conditions to obtain 
required test dynamic pressure, periodic calibrations of all the sensors 
were carried out using the portable shock tube apparatus, and the data 
were recorded on tape. Once this was performed, the side walls were 
closed and the model retracted into the quick insert chamber, in prepara- 
tion for the actual tests. When the proper reservoir conditions were 
achieved, the hydraulically actuated valve at the pebble bed heater was 
opened. As soon as steady test section conditions were reached, the 
model was injected into the stream for a short duration. During the 
model's stay in the stream for one minute, the data were taken on the 
20-channel tape recorder, both direct and FM modes of operation, at 
3.048 meters per second. Immediately after the run the recorded data were 
visually examined on an oscilloscope in order to ascertain that no data 
were lost. 


18 



At each wind-tunnel test condition four nearly identical tests were 
carried out in order to assure repeatability of tests and to increase the 
statistical confidence level of the recorded data during the analysis. 

At each test stream Mach number, visual flow observations were made 
to ascertain the location of the boundary-layer transition. The plate was 
sprayed with a saturated solution of Fluorene in a petroleum ether base 
before it was injected into the test section. Transition was found to 
occur 7 to 10 cm from the leading edge in all cases. This indicated that 
the inconel leading-edge portion and flat-plate surface did not match 
perfectly and thus acted as a triggering mechanism for the initiation of 
transition. This was an unexpected result, but assured a fully developed 
turbulent flow region over the instrximent plug. 


DATA ANALYSIS 

The analysis of the recorded wall pressure data was first examined 
for the stationarity conditions. The pressures recorded by sensors A to H 
are the fluctuations about a mean static pressure; and the root-mean- square 
(RMS) of these fluctuations from any one of the sensors for a given test 
stream condition remained an invariant irrespective of the length of record 
examined (above an optimum minimum playback of about 5 seconds at a tape 
speed of 3.048 m/sec) . The weak stationarity condition requires that first 
and second moments of time history of the pressure fluctuations be inde- 
pendent of time. These are satisfied by the results; namely, the time 
average of the fluctuating static pressure (about the mean wall static 
pressure) is zero and the RMS pressure is independent of time. 

The instrument noise was recorded with the entire model in the test 
section, with no stream flow in the wind tunnel and with the outputs of 
all the sensors recorded at the same amplifier gains as those in the actual 
Mach 7.4 tests. The noise results given in terms of ratios of an equi- 
valent RMS pressure for each sensor to dynamic pressure, q^, are in 
figure 20. The curve corresponding to each sensor is labeled. 
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MEAN AERODYNAMIC BOUNDARY-LAYER PARAMETERS 


, The aerodynamic parameters necessary for reducing the measured pressure 
fluctuations to generally accepted nondimens ional form are the boundary- 
layer thickness (<S), displacement thickness (6*), momentum thickness (0) 
and local skin-friction coefficient (Cf ) . In the tests in-situ measurement 
of these parameters was not possible due to malfunctioning of the boundary- 
layer traverse apparatus. However, Dies (ref. 45) who examined the existing 
literature on turbulent boundary-layer flows on a flat-plate model, derived 
an empirical formula in terms of Reynolds number, Re^ and Mach Nimiber, M^. 
These empirical relations are given below: 


X 


0.37 Re"® 

X 


^ 1 + 



(1.3 + 0.43 M^) (6/x) 

10.4 + 0.5 M^(l + 2xlO"®Re ) ° • ^ 3 3 

oo X 


( 1 ) 

( 2 ) 


= (1.3 + 0.43 M^) 


(3) 


Using these empirical relationships and the procedure given by Bies for 
obtaining local skin-friction coefficient, c^, the necessary boundary-layer 
parameters for present tasks were calculated from the known information on 
free- stream Mach number M^ and the Reynolds number at the center of the 
instrument plug, Re^^ and these are given in Table I. These results were ' 
compared with some available data on a similar flat-plate model tested in 
the Ames Research Center 3.5-Foot Hypersonic Wind Tunnel (ref. 46), and the 
agreement is good. In all the above considerations one has to keep in mind 
that our model is essentially a cold plate (i.e., surface temperature less 
than adiabatic wall temperature) in the test stream. 


RESULTS 

p Variation with 

The measured RMS surface pressures, p, from each sensor A to H are 
plotted against dynamic pressure, and are given in figures 21(a), (b) 
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(c) and (d) for the three Mach numbers; namely, 5.2, 7.4, and 10.4. 

Within the expected experimental scatter the RMS pressures obtained from 
the Kulite or the crystal sensor, located at the same downstream position 
from the leading edge, were found to be equal; that is, p^ = p^, Pg = Pp 
and so on. This indicates that the overall energy contribution from the 
high frequency portion of the power spectral density plots towards RMS 
pressure value is small. Hence, in all the discussions on the overall 
statistical properties of the pressures no distinction will be made 
regarding the two types of sensors. In figure 21 the RMS pressures are 
plotted against dynamic pressure, q^ and show that the pressures increase 
with the dynamic pressure, q^ and seem to reach a maximum around 
q^ = 5xlO‘*N/m^. A composite plot of all the RMS pressures shown in 
figure 22 is obtained from, figure 21(a) to (d) , in order to extract the 
general trend from the data. The RMS pressures from sensor A to D or 
from sensor E to H satisfies 


P 


A,E 




0.96 


0.92 : 0.89 


(4) 


and these are expressible in terms of the percent change of Reynolds 
number from sensor to sensor as compared with that of the leading sensor; 
namely, A or E. That is. 


^i " ^A,E ®^P(300R) 


(5) 


where 



( 6 ) 


i assumes A, B, C or D according to the sensor in consideration. Further, 
as Mach number increases, the value of p decreases at any given q^. 

The orientation of sensors with respect to mean stream direction is given 
in figure 5. In figure 22 the solid lines, obtained by the method of 
least squares, represent the measured data. An analytic expression has 
been obtained which fits all the measured data and is given below. 
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( 7 ) 


= [l.lxlO-"q^ - 6.959xlO-«q^ - 4.138]exp - [(l.OlTxlO"' 

+ 1.305xlO"®q - 1.261x10"^ ^q^)M 1 

This expression is valid for 5.2 < < 10.4 and for 6.5x10^ ^ ^oo £ 

7xl0‘*N/m^. The dashed lines were obtained from calculation of the analytic 
expression given in equation (7) for the sensor A or E. For convenience 
the ordinate scale at the right in figure 22 is given in db's based on 
reference RMS pressure equal to 2x10“ ®N/m^. 

p/p Variation with q^ 

In figure 23 (a) (b) , (c) and (d) the ratio of the RMS pressure to 

the mean static pressure at the surface, p/p is plotted against dynamic 
pressure, q^ for sensors A to H. The ratio of p/p decreases as q^ is 
increased from 2x10 ‘*N/m^ for M^ = 5.2 and 7.4 tests. But, for Mach 
number 10.4, this ratio of p/p, within the available data, increases as 
q^ is increased up to 2xio‘'N/m^. A composite plot obtained by separating 
the data for each Mach number 5.2, 7.4 and 10.4 is given in figure 24(a), 

(b) and (c) respectively. The composite plot (fig. 24) showsithat p/p 
reaches a maximum value around q = 2xlo‘*N/m^ for M =7.4 and then 
rapidly falls off to a small value as q^ is increased. Wind-tunnel 
operational limitations severely hampered efforts to obtain data for 
q^ < 2xlo‘*N/m^ at M^ = 5.2 or 7.4; and for q^ > 2xlO‘*N/m^ at 
M^ = 10.4. Thus, the entire range of q^ , for a given Mach number, was 
not possible. All the same it is of interest here to point out that the 
RMS pressures vary from a maximum of 25 percent to a minimum of 5 percent 
of the mean static pressure and depend on the dynamic pressure, q^. From 
the experimental results of Kistler and Chen (ref. 22), one notes the value 
of p to vary from 8 to 10 percent of the mean static pressures at 
M^ = 5, while Chen (ref, 21) reports for the same test conditions a vari- 
ation of up to 20 percent. Laufer (refs. 19 and 20) while examining the 
far field problem observed p/p equal to 1 percent outside the turbulent 
boundary layer, and pointed out the need for more experimental information 
in the near field at M^ > 5. Thus the present data seem to be more in 
agreement with Chen's data than with Kistler and Chen's results. 
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p/q Variation with q and M 

The ratios of the RMS pressure, p to the free-stream dynamic pressure, 
q^ are given in figures 25(a), (b) , (c) and (d) for each sensor and for 
all the three tested Mach numbers. These values are considerably smaller 
than the 6x10”^ value obtained for the subsonic , case (refs. 9 and 10). 
Further, the present data differ from the subsonic results in that the 
ratio P/'^ao depends on q^. Kistler and Chen (ref. 22) obtained P/^^ 
around 5x10“^ at M^ = 5 and the dependence of this ratio on was 

not evident from their data. In figure 26 we have reproduced the exact 
plots given in Chen's thesis (ref. 21, fig. 16). He gives P/^^ variation 
as a function of ^e^* and M^; and this is an order of magnitude larger 
than the data obtained from present tests (see fig. 28) . 

The composite plot of against q^ in figure 27 is obtained 

from replotting data from figure 25 from all the sensors. It is interesting 
to compare the data with that of figure 20, the instrument noise plot for 
each sensor and the associated electronics for M =7.4 tests. 

CO 

In order to assess the signal to noise in the instrumentation, we 

obtain the ratio of p/q^^ figure 27 and divide it by the converted 

noise value (from fig. 20 in equivalent P/q^ terms) at the selected q^ 

and M^. As an example consider sensor A, at M^ = 7.4 and q_^ = 5xlO‘*N/m^. 

In this case we obtain a signal to noise ratio of 2 . 5x10" Vl • 7x10“ '* = 14.7. 

Thus, the signal is nearly 15 times larger than the threshhold noise from 

all the electronic amplifiers in circuit for sensor A at q„ = 5x10 '* and 

M =7.4. 

00 

In figure 28, p/q„ is replotted against the free-stream Mach number, 

M for sensor A or E. The data shown inside the boxes are obtained from 
00 

figure 25(a). It can be seen as q^ increases the p/q^^ ratio tends 
towards Lowson's empirical curve (ref. 27). The data obtained in these 
tests generally fall between Houbolt (ref. 48) and Lowson's empirical 
curves and are smaller than the value obtained at subsonic speeds.* The 
data in the range of free-stream Mach number, 8.25 < M_^ < 9 are from 
Harvey et al, (ref. 37) who measured the velocity and temperature fluct- 
uations through the -use of electron-beam techniques. From these 


The vertical scale in these plots is logarithmic and does exaggerate the 
spread of data from the empirical curves. 
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measurements they were able to calculate the pressure fluctuations in- 
directly. These are shown in figure 28 and at high the data from the 

present investigation tends toward their results. In addition Harvey et al. 
calculated p/q^ values using Wallace's data for adiabatic and cold wall 
conditions and these are in figure 28 at = 9 . The adiabatic wall data 
of Harvey's coincides with that of Lowson's (ref. 47). The results from 
Chen (ref. 21) are an order of magnitude above all of other known experi- 
mental results. 

Variation of P/q„o a Function of and P/t^ as a Function of Reg 

In figure 29, p/q^ is plotted against the Reynolds niamber based on 
displacement thickness, Re^j*- The data satisfies the^ relation 

^ X 10^ = (40.587 - 1.678 M )Re"“-^“® (8) 

q “'6* 

^OO 


The variation of p/q^ as a function of suggests that one 

should examine these RMS pressures using local skin friction or local wall 

shear stress, as the characteristic nondimensionalizing parameter. 

Hence, we have plotted the RMS pressure data in terms of p/x against 

w 

the Reynolds number based on momentum thickness, Re». The data shown in 

w 

figure 30 may be represented by the relation 


= A Re® 2 
T 2 0 

W 


(9) 


where A 2 = 0.28 and = 0.27. 

The P/’’^^ data obtained in the present experiments for a given 

Re. ( =4.3x10'*) is presented in terms of M in figure 31. The data in 

general are lower in value than Kistler and Chen (ref. 22). Also, in 

figure 32 we have reproduced data from Chen's thesis (ref. 21). These are 

considerably higher than any other results. In general the present data 

show the same trends as Speaker and Ailman (re,f. 24) even though they are 

slightly higher in value. In the present investigation no special efforts 

were made to hold Re, constant in order to obtain p/x shown in fig- 

ure 31. But fortunately the test range covered a sufficient number of 

Re„ values from which the data for Re„ = 4.3x10'* could be selected for 
t) u 

all the three Mach numbers considered. 
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In figure 31 the data from the present investigation are extended to 
the subsonic range. This extrapolated value is in agreement with the 
results of Williams (ref. 31) and Belcher (ref. 32). 

Spectral Considerations 

The power spectrum of the recorded time history of the pressure fluc- 
tuations from each sensor A to H was carried out using a real time spectral 
analyzer and the associated electronic units described in the section on 
instrumentation. The entire spectral analyzer, a composite of several 
compatible electronic units, was calibrated by checking for the line spec- • 
tral output for sinusoidal inputs of known RMS voltage in the frequency 
range from 1 KHz to 1 MHz. The performance was also checked by examining 
the output of the analyzer for a square wave input at 1 KHz and 20 KHz 
and comparing the spectrum with the theoretical results. In addition, a 
random white noise generator of known RMS voltage with a frequency range 
up to 20 KHz was used as an input signal, and the spectral output obtained 
was compared with the known results. 

The nondimensional pressure spectral density and frequency were chosen 
to be of the form p^(f)U^/q^6* and f6*/U^ respectively. The choice of these 
parameters was due to the ready availability of the aerodynamic parameters 
and q^. Also, the physically meaningful boundary-layer thickness, 6* is 
retained. Each of the spectral curves represent 256 ensemble averages. 

The nondimensional pressure spectral density plots are given in fig- 
ures 33, 34 and 35 for the three Mach numbers investigated. The plots 
represent a composite of spectral data obtained by analyzing the data three 
times in frequency range of (i) 1 KHz £ f £ 20 KHz, (ii) 3.2 KHz £ f £ 

160 KHz, and (iii) 25 KHz £ f £ 320 KHz. In each of these plots there 
are eight curves, each representing the output of one particular sensor, 

A to H as indicated. On the first line at the top of each plot the test 
run number and the RMS pressures of sensors A to H (PRMS , N/m^) are given. 

The second line contains pertinent test condition information; namely, the 
free-stream dynamic pressure, q_^(QINF, N/m^), the stream velocity, 

(UINF, m/sec), Reynolds number in millions, Re^ (Reynolds number), and the 
displacement thickness 6* (DELSTR, m) . Using the Reynolds number Re^ 
given in the plots and the aerodynamic parameters 6, 6*, 6 and c^ 
given in Table I, any of the commonly accepted nondimens ionalizing parameters 
can be obtained for replotting the presented data. 
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The RMS pressures in figures 33-35 were obtained using Simpson's rule 
for area integration of the spectral density plot and were checked against 
the values that were directly obtained through the use of a Ballantine 
meter. The agreement between these two ways of obtaining p is within 
2 percent of the reduced pressure data. 

One common noticeable feature of the spectral density plots is that a 
sharp drop occurs at large Strouhal number, f6*/U^. This characteristic 
drop is attributable to the finite size of the sensor. As the turbulent 
eddy size becomes comparable or smaller than the sensor diameter, the 
sensor response falls off drastically. In a convected turbulent boundary 
layer flow the small eddies correspond to high frequency components. Thus, 
at some high frequency (determined by the ratio of sensor diameter to 
boundary layer displacement thickness, d/6*) cancellation between positive 
and negative pressure fluctuations begin to occur over the sensor surface. 
Willmarth (ref. 12), Serafini (ref. 18) and others have suggested on the 
basis of their experimental investigations that the ratio of d/6* is to 
be kept as small as possible in order to obtain high frequency (small eddy 
size) resolution of the convected pressure field. In figure 36 the effect 
of the ratio of sensor diameter d , to displacement thickness, 6* is 
indicated. In addition, as the size of the sensor is reduced the sensiti- 
vity is decreased. In spite of all these considerations, the finite physical 
size of the sensor allows us to record only the pressure that has suffered 
some attenuation through the filtering effect of the measuring instrument. 

In order to obtain the true value from the measured quantity a correction 
is essential if d/6* is large. Corcos (ref. 49), Willmarth and Roos 
(ref. 14), Chandiramani (ref. 50), Gilchrist et al . (ref. 51) and White 
(ref. 52) investigated this problem both experimentally and theoretically. 
From their results one is able to obtain the expected reduction in pressure 
sensitivity of the sensor due to its finite size; and this is given in fig- 
ure 37 . In the present investigation using an average free-stream velocity 
of 1.4 Km/sec the sensors satisfy the following condition. 


Kulite Sensor : 

0.10 

A 

A 
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From figure 12, the calculated attenuation for the crystal sensor used in 
the tests at 300 KHz is equal to 0.3 dbs. Thus, the corrections for the 
size of the sensor are negligible and hence no corrections were made. 

For a free-stream Mach number of 5.2, the power spectral density 
exhibits some scatter at low Strouhal number between adjacent measuring 
stations. Similar results have been observed by previous investigators 
while carrying out tests in the supersonic region of = 3.5 (refs. 24 

and 27). However, at higher Strouhal numbers all these data coalesce 
nicely. Various explanations have been offered for low Strouhal number 
behavior; namely, some weak shock wave interaction, random eruptions of 
turbulence bursts due to surface roughness and the probability of having 
a region where the boundary layer turbulence is not fully developed. In 
the present tests one finds that the pebbles in the heater chamber are 
literally lifted up and settled back for the high dynamic pressure and mass 
flows at a test section Mach niamber of 5.2. This process releases all the 
disintegrated alumina oxide particles as dust. In order to alleviate the 
problem a periodic wind-tunnel dust blow-down, without the model in the 
stream was carried out. In spite of this precaution the problem of dust 
was present for the Mach 5.2 tests. 

The phenomenon of dust in the stream disappears at Mach numbers of 
7.4 and 10.4 and is one of the unexplained idiosyncrasies of the tunnel 
operation. The effect on the power spectral density plot is less scatter 
of data at these Mach numbers (M = 7.4 or 10.4) . 

CO 

In order to obtain the general trends and to arrive at an analytic 
expression, a composite of several of these power spectral density plots 
for a given Mach number were examined. The spectral density is in general 
satisfied by the following equation for all the Mach numbers considered. 
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The constant A is related to the RMS pressure as can be seen by simple 
integration. 
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An analytic expression for which fits most of the measurements has 

already been obtained (eq. (7)). Thus, all the power spectral density 
plots for the pressure fluctuations obtained in these tests can be collapsed 
by considering 


i,e. 
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Lowson (ref. 47) suggested a similar expression after examining some of 
the supersonic experimental results. The amplitudes of the power spectral 
density from the plots presented are in general lower than the measurements 
given by Chen (ref. 21) Kistler and Chen (ref. 22), Coe (refs. 27 and 28) 
and others. The data presented by Speaker and Ailman (ref. 24) for 
attached turbulent boundary layer flow agrees with results of this investi- 
gation. This can be predicted from p/q^ plots versus Mach numbers (see 
fig. 32) which indicates lower p/'l.^o values as M^ is increased for a 
given q^ value. 


Correlation Functions 

The correlation functions of interest in the present work are the 
space-time and auto correlations of the recorded surface pressure fluctu- 
ations. The space-time correlation is defined as 
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where the correlation function is normalized by the use of appropriate RMS 
pressures. Two-point space-time correlations provide information on the 
convective velocity of the pressure fluctuations and on the decay rate of 
the pressure signature as it is convected downstream. 

The auto-correlation function is a special case of the space-time 
correlation and is obtained by letting the space separation vector ^ 
between sensors go to zero. This essentially implies that the pressure 
fluctuation considered is the output of one sensor with the usual time 
displacements carried out for obtaining the correlation function. In this 
case, since 5=0, i = j. For a stationary process the auto correlation 
can also be obtained by considering the Fourier transform of the power 
spectral density function.* This can be expressed by using equation (10) 
and considering its Fourier transform. Hence 




eiBf*df* 


(1 + f *^) ^2 


where 


U t 

f* = , B = 2tt 

U ■ ° 


and q A = p 


and finally 


R(t) 


00 

-I 


iBf^ 


(1 + f*^) V2 


df * 


(17) 


(18) 


(19) 


This integration can be carried out by contour integration methods. 

The entire correlation analysis has been carried out for the free- 
stream Mach number of 7.4 using the correlator described in the section 


*The pressure fluctuation data satisfy the stationarity condition;, and thus 
the auto correlation the power spectral density function, 

p^(f) are Fourier transforms of each other. 



on instrumentation. This analysis is not as exhaustive as the conside- 
rations of RMS pressures and the power spectral densities. In addition, 
in this analysis only broad frequency bands were considered (1 KHz £ f 
£ 300 KHz) . The pressure sensors are located a fixed distance apart 
in the stream direction in order, to measure the longitudinal statistical 
characteristics of pressure fluctuations. No lateral correlation consi- 
deration was undertaken in this investigation. 

Two space-time correlations can be analyzed at one time with the 
correlator, which accepts three separate pressure data as input simul- 
taneously and yields output voltages proportional to the correlation 
function, for instance, order 

to obtain all the three space-time correlations of interest with this 
correlator the analysis had to be performed twice for each set of test 
conditions. Care was exercised to examine the data in the same time 
interval for both the correlation analyses. Thus the correlations ,x) 

and R (C ,t) from analyzing the data the first time, and R ,x) and 

O A^ O 

R (C ,t) from the second analysis were obtained. In figure 38 only one 
AL/ O 

trace of the common R (5 ,t) correlation is given. 

Av-' O 

The convective velocity can be defined in two ways, each giving a 
slightly different result. The, two convective velocities are defined as 



where x corresponds to a x that satisfies the condition 


3R(g,T) 

3x 


( 21 ) 


and 



where x corresponds to a x that satisfies the condition 


3R(^,x) 

3|C| 


0 • 


( 22 ) 


(23) 
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From the present analysis only can be obtained since the sensors 

fixed in the instrument plug cannot be easily varied in their separation 
distances from each other. Thus one measures the time delay, t, required 
between the peak values of two consecutive correlations , for instance 
and or and and since the displacement distance between sensors 

are known we obtain the so-called convective -velocity as equal to 5 /t • . 

The convective velocities so obtained are in the range 

0.60 £ -yi < 0.70 
00 

and these values agree reasonably well with measurements given by Kistler . 

and Chen (ref. 22), Chyu and Hanly (for Cj/<S* = 1.5 - 2.0,- ref. 53) and 

Bull (ref. 15) . Even though a systematic space-time correlation for the 

other two Mach numbers = 5.2 arid 10.4 were. not performed, the results 

indicate Kistler and Chen's (ref. 22) findings of the /U^ ratio as a.' 

function of 'to be valid; namely, 0^ /U^ decreases from a value of 

0.8 at subsonic flow to 0.65 at supersonic and remains constant beyond 

M = 3.5. Similar results have been obtained by Coeri(ref. 27).. From an 
00 " . 

examination of shadowgraph pictures Jedlick et al. (ref. 54) reported con- 
vective velocities of. bursts of turbulence in the boundary layer to be 
about 60 percent of the free-stream velocity. 

From the space-time correlation considerations Uq /U^ = 0.65 was 
obtained for a separation distance of = 0.95 cm, and Uq /U^ = 0.700 

for a separation distance of 2?^ = 1.9 cm. Thus, as the separation 
distance increases, there is a slight increase in Uc This is in 

agreement with the work of Bull (ref. 15) who reports an increase in the 
/U^ value as separation distance increases for the subsonic case. 

These results indicate that the small scale (or large wave number) 
components of pressure fluctuations travel downstream slowly and lose 
their identity rapidly; while the large scale (or small wave number) 
components are convected rapidly and lose their identity slowly. 

The auto-correlation plot. is useful if one is interested in the micro- 
scale of pressure fluctuations (analogous to a velocity field) . A repre- 
sentative auto-correlation plot is given in figure 39. The vertical scale 
in these plots are relative and are the voltage outputs of the correlator. 

An osculating parabola is drawn in dashed lines in order to obtain the 
microscale length of the pressure fluctuating field. The average of the 
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intercept of this curve on the abscissa is around 50 psec. If one uses 
the average convective velocity of 900 m/sec then the microscale is equal 
to 50x 10~ ^ X 900x 100 = 4.5 cm, and is greater than the distance between 
sensors . 

CONCLUSIONS 

Fluctuating pressures have been measured on a flat-plate model at 
hypersonic Mach numbers of 5.2, 7.4 and 10.4 beneath an attached turbulent 
boundary layer for a range of Reynolds numbers from 1 to 15 millions and 
for a dynamic pressure range of 6800 to 66000 N/m^ . The following con- 
clusions can be drawn from the analysis. 

1. With proper care the piezoelectric transducers could be success- 
fully used in a pebble heated hypersonic wind tunnel, for measuring surface 
pressure fluctuations. 

2. Over the range of Mach number tested, the RMS pressures increase 
with the dynamic pressure and reach a plateau. For any dynamic pressure, 
the RMS pressure decreases with an increase in Mach number. The pressure 
intensity decreases systematically as the pressure signature is convected 
downstream from sensor to sensor. 

3. The variation of the RMS pressures normalized with mean static 
pressure is 8 to 25 percent, and the peak values occur at. dynamic pressures 
around 2xlO'^N/m^. 

4 . The ratios of RMS pressures to dynamic pressure all fall below 
4x10“^ for these hypersonic wind-tunnel tests and this ratio decreases 

as dynamic pressure or Mach number is increased. For high dynamic pressure 
level at hypersonic speeds, the normalized RMS pressure to dynamic pressure 
could be approximated by Lowson's empirical formula (ref. 27). 

5. The ratio of RMS pressure to dynamic pressure varies as (Re^*‘^) 
for all three test Mach numbers. 

6 . The ratio of RMS pressure to wall shear stress falls between 2 
and 3 and considerably lower than Kistler and Chen's results (ref. 22). 

As Mach number increases, this ratio decreases in accordance with the 
trends observed by Speaker and Ailman (ref. 24). 
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plots 


7. The nondimensionalized power spectral density p^(f)U^/q^6* 
against frequency parameter f6*/U^ show that all the plots can be 
collapsed to a single analytic representation given in equation (15). 

The contribution from high frequency components (100-300 KHz) of the 
pressure fluctuations to the pressure intensity (p) is very small. Even 
though the amplitudes of the high frequency components are small , these 
might be of interest for structural response considerations leading to 
fatigue problems. 

8. The space-time correlation considerations yield convective 
velocities that are equal to 0.7 of free-stream velocities. 

Nielsen Engineering & Research, Inc. 

Mountain View, California 
June 1973 
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18.0 , 

128 

7.4 

1.487 

2.413 

1,855 

1132.2 

3.723 

1.34 

14.53 

9.78 

0.396 

18,0 

129 

7.4 

1.460 

2.413 

1.855 

1098.9 

3.585 

1.34 

14.53 

9.78 

0.396 

18.0 

130 

7.4 

1.460 

2.413 

1.855 

1098.9 

3.585 

1.34 

14.53 

9.78 

0,396 

18.0 

122 

7.4 

1.423 

1.724 

1.310 

1040.0 

2.965 

1.13 

15.16 

9.91 

0,409 

19.8 

123 

7.4 

1.463 

1.724 

1,296 

1101.1 

3.172 

1.13 

15.16 

9.91 

0.409 

19.8 

124 

7.4 

1.451 

1.724 

1.255 

1084.4 

3.172 

1.13 

15.16 

9.91 

0.409 

19.8 

125 

7.4 

1.489 

1.724 

1.241 

1135.0 

2.896 

0.93 

17.47 

10.51 

0.423 

22.0 

110 

10.4 

1.402 

1.448 

9.722 

989.4 

2.344 

4.18 

11.84 

8.79 

0.185 

14.2 

111 

10.4 

1.463 

1.379 

9.625 

1066.7 

3.309 

3.47 

11.73 

9.24 

0.195 

15.2 

112 

10.4 

1.512 

1.379 

9.625 

1122.2 

2.689 

3.14 

11.99 

9.40 

0.198 

16,0 

113 

10.4 

1.512 

1.379 

9.625 

1122,2 

2.689 

3.14 

11.99 

9.40 

0.198 

16.0 

114 

10.4 

1.494 

1.379 

9.687 

1103.3 

2.482 

3.25 

11.89 

8.99 

0.195 

15.8 

114 

10.4 

1.524 

1.379 

9.687 

1142.2 

2.404 

3.14 

11.99 

9.35 

0.198 

16.0 

115 

10.4 

1.530 

1.379 

9.687 

1152.2 

1.861 

3.03 

12.09 

9.40 

0.199 

16.2 

116 - 

10.4 

1.402 

0.758 

4.227 

986.1 

1.793 

1.82 

13.72 

10.54 

0.221 

20.1 

117 

10.4 

1.497 

0.689 

4.220 

1115.0 

2.965 

1.43 

14.53 

12.14 

0.231 

22.7 

118 

10.4 

1.494 

0.689 

4.227 

1108,9 

2.413 

1.43 

14.53 

12.14 

0,231 

22.7 

119 

10.4 

1.494 

0.689 

4.220 

1107,2 

1.930 

1.43 

14.53 

12.14 

0.231 

22.7 

120 

10.4 

1.527 

0.689 

4.213 

1192.8 

2.413 

1.23 

15.14 

12.45 

0.238 

24.0 


Table I.— Summary of test conditions ond oerodynomlc parameters. 
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NOZZLES 


PEBBLE BED HEATER 


(a) Overall view of the test facility 


VACUUM SPHERES 


HIGH PRESSURE 
STORAGE 


COOLING TOWERS 


ARER COOLER 


TEST SECTION 
CONTROL 


AUX. EQUIPMENT 
BUILDING 


CRYOGENIC SEPARATION 
TOWER 



(b) Pebble bed heater, nozzle and test section 


Figure I Ames 3.5-Foot Hypersonic Wind Tunnel. 



o 



Figure 2.- Model injector, test section, and 
interchangeable nozzles. 


•J * 








M 



Figure 4.- Test model and protective enclosure 
for the instruments. 




O'rings location for 

SHOCK MOUNTING 


Kel-F thermal 
INSULATION 


A-D - Kulite Sensor 

E-F - Crystal Sensor 

- Thermocouples 

LI I - Static Pressure 

Orifice 

tj, - 9.525 nmi. 


Figure 5 - Instrument plug showing various sensors 



iHERMOCOUPLE 


RING FOR SOURCE' 


FOLLOWER 


Leads from Kulite sensor 


Crystal se nsor 


Thermocouple 


From static pressure 

OR I P I C & 


SoURCE-rOLLOWLk 


Figure 6.- Underside of the instrument plug showing sensors 
and signal conditioning amplifiers. 
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Figure 7.- Shock-tube calibration apparatus. 


0.5 mm 

CRYSTAL SENSOR 


RTV INSULATION 


3 mm diam 



Figure 8.- Schematic diagram of crystal sensor. 
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3 


Sensitivity, iJv/(N/m^) 

sitivity of all pressure sensors in decibel units. 




Ln 

o 


BULLET 
M ^ 3.0 



CRYSTAL 

SENSOR 


Figure 11 .-Determination of time rise of crystal sensor. 
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Frequency, Hz 

Figure 12.- Crystal sensor characteristics - amplitude and phase versus frequency 





6.5x10''- 1 lx 10'' N /m 


550 - 1400° K 


44 cm 


INSTRUMENTED 
— PLUG ( 
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V^A 




QUICK INSERT 
MECHANISM 


MODEL LOCATION / 
PRIOR TO INSERTION 




Figure 13.- Flat-plate model installed in Ames , 3 . 5-Foot Hypersonic Wind Tunnel 





Surface Temperature (°C) 
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Figure 1*^.- Surface temperatures on flat-plate model during tests 




To the tape 
recorder 
channel 


0.22 pf condenser for blocking DC voltage. 

10 leakage path resistance for charge 
build-up in condenser. 

PAR CR-4 amplifier operation-single ended. 


(a) Schematic wiring diagram for Kulite sensor for 
unsteady pressure measurement. 


Model 402A11 
Piezotronic 
Source- 
Follower 



To the tape 
recorder 
channel 


0.22 yf - DC blockage capacitance. 

MCL 1302 - 2 ma constant current diode. 

2500n - resistance used in order to obtain less than 
1 db noise figure from CR-4 amplifier. 

(b) Schematic wiring diagram for crystal sensor. 


Figure 15.- Electrical circuit used for pressure sensors. 
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SOURCE IMPEDANCE AT 300°K-OHMS 






Relative Voltage Gain 


Relative Voltage Gain = Vo°lta?e -) (cilir) 1 


Astrodata Model 886 
(8-amplifier data) 

Princeton Applied Research PAR CR-4 
(10-amplifier data) 




/ '/ Roll-off set 
at 10 Hz 


-6 dB Point, 


\ Roll-off 
set at 
300 KHz 


Roll-off set at 300 Hz 


Frequency, Hz. 

Figure 17-- Amplifier frequency response characteristics 





Lead Wire (First 
Thermocouple Element) 



Form Thermoelectric 
Junction 


Figure 19.- Cross-sectional schematic of 
TCS Series Probe. 









Free-stream dynamic 
(c) Pressure data 
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Figure 21 .- 


Root-raean-square of surface pressure fluctuations, p(N/m 
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Root-mean-square of surface pressure fluctuations, p(N/m 



(dB's, reference pressure 2x10“ ^N/m 
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Ratio of RMS pressure to dynamic pressure, p/q 
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Figure 28.- 


468 
Free- stream Mach number, M 


Ratio of RMS pressure to dynamic pressure 
plotted against Mach number. 




10 “* 


2x10“* 4x10“* 6x10“* 8x10“* 10^ 

Reynolds number based on displacement thickness, 

Figure 29.- Normalized FIMS pressure, p/q^ against Reynolds number, 

-a 
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Figure 31.- Ratio of RMS pressu 

against Mach 
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(a) q = 50332 N/m^ and Re = 4.55x10®. 

\ • *oo 

Figure 33.- Power spectrum of pressure fluctuations, = 5.2. 
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pMf)U 



(c) = 35163 N/m^ and Re^ = 2.1x10®. 

Figure 33.- Concluded. 
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(f)u 


abcdefoh 

RUN NO. 107-I6-00-C0 FRMS(N/M2) = 126.86 114.45 122.05 108.24 126.86 114.45 122.05 108 24 

1 OKlO"**^ alNF(N/M2) : 55158. UINFCM/S) : 1254.4 REV. MO. =8.9 9 (MILS) DELSTR(M) = .0064 



f6*/U„ 

(b) = 55158 N/m^ and Re^ = 8.99x10®. 

Figure 34.- Continued. 
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P^(f)U„ 



(c) q = 46884 N/m^ and Re = 4.33x10®. 

«> X 

Figure 34.- Continued. 
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(f)u 


ABC0EF6H 
RUN NO. 115-lS-CO-CC PRMS(N/M2) = 40.6T 39. SO 34.47 31.71 40.67 39.30 34.47 31.71 

1 OXlO*^ 9INF<N/M2l = 137B9. UlNF<N/5) r 1330.1 RE7. NO. : 3l0 3 (NILS) OELS7R(NI z .0093 



(a) = 13789 N/m^ and Re^ = 3.03x10®. 

Figure 35.- Power spectrum of pressure fluctuations, = 10.4. 
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Time delay, t (y sec) 
Figure 38.- Space-time correlation for 
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